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Globally structured 3D analysis-suitable 
T-splines: definition, linear independence and 
m-graded iocai refinement 


Philipp Morgenstern* 
May 12, 2016 


This paper addresses the linear independence of T-splines that correspond to refinements 
of three-dimensional tensor-product meshes. We give an abstract definition of analysis- 
suitability, and prove that it is equivalent to dual-compatibility, wich guarantees linear 
independence of the T-spline blending functions. In addition, we present a local refinement 
algorithm that generates analysis-suitable meshes and has linear computational complexity 
in terms of the number of marked and generated mesh elements. 

Keywords: Isogeometric Analysis, trivariate T-Splines, Analysis-Suitability, Dual-Compatibility, Adap¬ 
tive mesh refinement. 


1. Introduction 

T-splines lUl have been introduced as a free-form geometric technology and were the first tool of interest 
in Adaptive Isogeometric Analysis (IGA). Although they are still among the most common techniques 
in Computer Aided Design, T-splines provide algorithmic difficulties that have motivated a wide range 
of alternative approaches to mesh-adaptive splines, such as hierarchical B-splines l|2j|3l, THB-splines 
|4l, LR splines ||5l, hierarchical T-splines ||6l, amongst many others. 

One major difficulty using T-splines for analysis has been pointed out by Buffa, Cho and Sangalli 
fT*! , who showed that general T-spline meshes can induce linear dependent T-spline blending functions. 
This prohibits the use of T-splines as a basis for analytical purposes such as solving a discretized partial 
differential equation. This insight motivated the research on T-meshes that guarantee the linear inde¬ 
pendence of the corresponding T-spline blending functions, referred to as analysis-suitable T-meshes. 
Analysis-suitability has been characterized in terms of topological mesh properties p8i| and, in an alter¬ 
native approach, through the equivalent concept of Dual-Compatibility f9i|. While Dual-Compatibility 
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has been characterized in arbitrary dimensions ||T(1 . Analysis-Suitability as a topological criterion for 
linear independence of the T-spline functions is only available in the two-dimensional setting. 

In this paper, we introduce analysis-suitable T-splines for those 3D meshes that are refinements of 
tensor-product meshes, and propose an algorithm for their local refinement, based on our preliminary 
work in ifTTl . In addition, we generalize the algorithm from ifTTl by introducing a grading parameter 
m that represents the number of children in a single elements’ refinement. This allows the user to fully 
control how local the refinement shall be. Choosing m large yields meshes with very local refinement, 
while a small m will cause more wide-spreaded refinement. The former yields a smaller number of 
degrees of freedom, while the latter reduces the overlap of the basis functions and hence provides 
sparser Galerkin and collocation matrices. 

This paper is organized as follows. Section [^defines the initial mesh and basic refinement steps and 
introduces our new refinement algorithm. Sectionj^then characterizes the class of ‘admissible meshes’ 
generated by this algorithm. In Section we give a brief definition of trivariate odd-degree T-splines. 
In Section]^ we give an abstract definition of Analysis-Suitability in the 3D setting and prove that all 
admissible meshes are analysis-suitable. In Sectionj^we define dual-compatible meshes, and prove that 
analysis-suitability and dual-compatibility are equivalent, and that all dual-compatible meshes provide 
linear independent T-spline functions. (Figure [T] illustrates this “long way” to linear independence.) 
Section proves linear complexity of the refinement procedure, and conclusions and an outlook to 
future work are finally given in Section 
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Figure 1: How we prove linear independence of the T-splines induced by the generated meshes. 


2. Adaptive mesh refinement 

This section defines the new refinement algorithm and characterizes the class of meshes which are gen¬ 
erated by this algorithm. The algorithm is essentially a 3D version of the one introduced in lUTlI . with 
the additional feature of variable grading. The initial mesh is assumed to have a very simple structure. 
In the context of IGA, the partitioned rectangular domain is referred to as index domain. This is, we 
assume that the physical domain (on which, e.g., a PDF is to be solved) is obtained by a continuous 
map from the active region (cf. Section |^, which is a subset of the index domain. Throughout this 
paper, we focus on the mesh refinement only, and therefore we will only consider the index domain. 
For the parametrization and refinement of the T-spline blending functions, we refer to IIT2ll . 
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Definition 2.1 (Initial mesh, element). Given X,Y,Z e N, the initial mesh Qq is a tensor product mesh 
consisting of closed cubes (also denoted elements) with side length 1, i.e., 

Qo := jU-l,^]xb;-l,3;]x[z-l,z] |xe{l,...,l),je{l,...,F},z€{l,...,Z)). 

The domain partitioned by Qq is denoted by Q. := (0, X) x (0, Y) x (0, Z). 

The key property of the refinement algorithm will be that refinement of an element K is allowed 
only if elements in a certain neighbourhood are sufficiently fine. The size of fhis neighbourhood, 
which is denofed (p, m)-pafch and defined fhrough fhe definitions below, depends on fhe size of K, fhe 
polynomial degree p = {pi,P 2 ,P 3 ) of fhe T-spline functions, and fhe grading paramefer m. For fhe 
sake of legibilify, we assume fhaf pi,p 2 ,P 3 are odd and greafer or equal fo 3. (For commenfs on even 
polynomial degrees, see Secfionj^) 

Definition 2.2 (Level). The level of an element K is defined by 


{(K) :=-log,J^|, 


where m is the manually chosen grading parameter, i.e., the number of children in a single elements’ 
refinement, and |.K'| denotes the volume of K. This implies that all elements of the initial mesh have 
level zero and that the refinement of an element K yields m elements of level £{K) + 1. 

Definition 2.3 (Vector-valued distance). Given x e Q. and an element K, we define their distance as the 
componentwise absolute value of the difference between x and the midpoint of K, 

Dist(/r, x) := abs(mid(/f) - x) e 
with abs(y) := (|yi|,|y2l,|y3l)- 

For two elements K\,K 2 , we define the shorthand notation 


Dist(A'i,.K' 2 ) := abs(mid(A'i) - mid(A' 2 )). 


Definition 2.4. Given an element K, a grading parameter m > 2 and the polynomial degree p = 
{P\,P 2 , P 3 ), we define the open environment 


where 


U^’'"{K) :={x€n\ Dist(/:,;c) < 


DP’'"(k) := 


m 


(pi + hP2 + hp3 + 5 ) 




m 


m 5 ^ ■ 2 

-(k-2)/3 (' Pi+3/2 P2+3/2 


,P3 + i) 


if k = 0 mod 3, 
if k = 1 mod 3, 
if k = 2 mod 3. 


The (p, m)-patch of K is defined as the set of all elements that intersect with environment of K, 


:= {K' eg\K' D U^’'^{K) + 0). 

Note as a technical detail that this definition does not require that Keg. See also Figure for 
examples. 
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Figure 2: Examples for the (p, m)-patch of an element K, for p = (3,3,3), m = 3 and €{K) = 2,3,4. 


Remark. By definition, the size of the (p, m)-patch of an element K scales linearly with the size of K 
and with the polynomial degree p. Since DP’'”(^) is decreasing in m, choosing m large will cause small 
(p, m)-patches and hence more localized refinement. 

In the subsequent definitions, we will give a detailed description of the elementary subdivision steps 
and then present the new refinement algorithm. 

Definition 2.5 (Subdivision of an element). Given an arbitrary element K = [x, x+x] x [y, y+y] x [z, z+2], 
where x, y, z, x, y, z € K. and x, y, z > 0, we define the operators 

subdivx(^) := I [x+ ^x,x+ ^x] x [y,y+ y] x [z,z + z] I j € {l,...,m)), 

subdiVyC.^) := I [x,x + x] X [y + ^y.3' + x [z,z + z] | y e {1,... ,m)), 

and subdiv^C/sT) := | [x, x + x] x \y,y + y] x [z + ^z,z + ^z] | y € {1,..., m}}. 


These operators will be used for x-, y-, and z-orthogonal subdivisions in the refinement procedure. 
Their output is illustrated in Figure 


Definition 2.6 (Subdivision). Given a mesh ^ and an element K e Q, we denote by subdiv(@^, K) the 
mesh that results from a level-dependent subdivision of K, 


subdiv(@, K):=g\ {/«:) U child(/:), 


subdivx(.^0 


with 


child(^) := - 


subdivy(.^0 


subdivz(.^0 


if W 
if W 
if W 


0 mod 3, 

1 mod 3, 

2 mod 3. 


Definition 2.7 (Multiple subdivisions). We introduce the shorthand notation subdiv(^, Af) for the sub¬ 
division of several elements Af = {Ai, .. .,Kj} c defined by successive subdivisions in an arbitrary 
order, 

subdiv(^, Af) := subdiv(subdiv(... subdiv(^, j). 

We will now define fhe new refinement algorithm through the subdivision of a superset closP’'”(^, Af) 
of the marked elements Af. In the remaining part of this section, we characterize the class of meshes 
generated by this refinement algorithm. 
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Figure 3: Elementary subdivision routines for m = 3: ;c-orthogonal subdivision of an element with 
level 0 (left), y-orthogonal subdivision of an element with level 1 (middle), and z-orthogonal 
subdivision of an element with level 2 (right). 


Algorithm 2.8 (Closure). Given a mesh Q and a set of marked elements A1 c ^ to be refined, the 
closure clos^’'”{0, Al) of A1 is computed as follows. 

M:=M 

repeat 

for all A e Al do 

M--M^[K'€ I {{K') < {(K)} 

end for 

until Al stops growing 
return c1osP’'"(@, Al) = Al 

Algorithm 2.9 (Refinement). Given a mesh g and a set of marked elements Al £ ^ to be refined, 
refP’™!^, Al) is defined by 

reP’'”(g,M) := subdiv(@,closP’'”(@, Al)). 

An example of fhis algorifhm is given in Figure]^ 

Example 2.10. Consider an inifial mesh fhaf consisfs of 4 x 5 x 8 cubes of size 1x1x1. We refine fhe 
mesh by marking fhe lower lefl fronl comer elemenf repeafedly until if is of fhe size x x The 
resulfing meshes for differenf choices of m are illusfrafed in Figure]^ and fhe resulfs are lisfed below. 


Figure 

m 

number of 
refinement steps 

number of 

new elements 


5a 


2 

12 

10728 


5b 


4 

6 

3175 


5c 


16 

3 

1030 


3. Admissible meshes 

In fhe subsequenf definilions, we infroduce a class of admissible meshes. We will fhen prove fhaf fhis 
class coindices wifh fhe meshes generated by Algorifhm |2.9| 

Definition 3.1 ((p, m)-admissible subdivisions). Given a mesh g and an element Keg, the subdivision 
of K is called {p,m)-admissible if all K' e satisfy £{K') > €{K). 

In the case of several elements At = {Ai,...,A/} c g, the subdivision subdiv(^,At) is (p,m)- 
admissible if there is an ordering (cr(l),..., cr{J)) (this is, if there is a permutation cr of {1,..., 7)) such 
that 

subdiv(^, At) = subdiv(subdiv(... subdiv(^, Ao-(i)),...), Ao-(y)) 
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Figure 4: Example for Algorithm |2.9[ with p = (3,3, 3),m - 3 and A1 = {A') with £{K) = 2. In the first 
iteration of the for-loop, all coarser (level 1) elements in the (p, m)-patch of K are marked 
as well, in the second iteration, all coarser (level 0) “neighbours” of those elements are also 
marked. Since there are no elements that are coarser than level 0, the third iteration does not 
change anything. Hence the for-loop ends, and all marked elements are subdivided in the 
directions that correspond to their levels. 


is a concatenation of (p, m)-admissible subdivisions. 

Definition 3.2 (Admissible mesh). A refinement 0 of Qq is (p, m)-admissible if there is a sequence of 
meshes = Q and markings Mj Q Qjfov j - 0,... ,J - such that = subdiv(^y, Al/) 

is an (p, m)-admissible subdivision for all j = 0,..., 7 - 1. The set of all (p, m)-admissible meshes, 
which is the initial mesh and its (p, m)-admissible refinements, is denoted by For the sake of 

legibility, we write ‘admissible’ instead of ‘(p, n 2 )-admissible’ throughout the rest of this paper. 

Theorem 3.3. Any admissible mesh Q and any set of marked elements satisfy 


refP’'”(@,M) € AP’'”. 


The proof of Theorem 3.3 given at the end of this section relies on the subsequent results. 


Lemma 3.4. Given an admissible mesh Q and two nested elements K Q K with A, A € IJ AP’'”, the 
corresponding (p, m)-patches are nested in the sense c ^P’'”(A). 


The proof is given in Appendix |A.1| 

Lemma 3.5 (local quasi-uniformity). Given K e Q e AP’"*, any K' € ^P’'”(A) satisfies i{K') > 7(A)-1. 
The proof is given in Appendix |A.2| 

Proof of Theorem [!0] Given the mesh Q € AP’'” and marked elements Af £ ^ to be refined, we have 
to show that there is a sequence of meshes that are subsequent admissible refinements, with Q being 
the first and refP’^f^, At) the last mesh in that sequence. 
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(a) m - 2 


(b) m — 4 


(c) m — 16 


Figure 5: Refinement examples for p = (3,3,3) and different choices of m. In all cases, the initial mesh 
consists of 4 X 5 X 8 cubes of size 1x1x1, and is refined by marking the lower left front 
comer element repeatedly until it is of the size x x ■ 


SetAl := closP’'”(@,M)and 


L:=max^(Af), L:=min^(Af) 

Mj ■-[K^M\ i{K) = j] for j = L,...,L 

9l := Q, 9j+i '■= subdiv(6^y, Mj) for j ^ L,..., L. (1) 


It follows that refP’'”(^, At) = 9z+i- "'ll! show by induction over j that all subdivisions in ([T]) are 

admissible. 

For the first step j = L, we know {K' € Af | i{K') < L) = 0, and by construction of Af that for 
each K e AIl holds {K' e Q'^'’”{K) \ €{K') < i{K)} c Af. Together with i{K) = L, it follows for any 
K € Ml that there is no K' € with i{K') < €{K). This is, the subdivisions of all K € Af^ are 


admissible independently of their order and hence subdivf^i, AIl) is admissible. 

Consider an arbitrary step y e {L,..., L) and assume that Ql, ■ ■ ■ ,9j are admissible meshes. Assume 
for contradiction that there is A € Af/ of which the subdivision is not admissible, i.e., there exists 


K' € QY'\K) with €{K') < i{K) and consequently K' ^ Af, because K' has not been refined yet. It 
follows from the closure Algorithm 2.8 that K' ^ Q. Hence, there is K eQ such that K' c K. We have 
{{k) < i{K') < i{K), which implies i{k) < £{K) - 1. Note that K e Q because Afy Q M Q Q. From 
K' € g^f'^iK), it follows by definition that K' n f/P’^fA) 0, and K' c K yields K n f/P’”XA) t 0 and 
hence K € ^P’^fA). Together with {(K) < i{K) - 1, Lemma 3.5 implies that Q is not admissible, which 
contradicts the assumption. □ 
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4. T-spline definition 


In this section, we define trivariate T-spline functions corresponding to a given admissible mesh. We 
roughly follow the definitions from ifTTlI . 

Definition 4.1 (Active nodes). For each element K - {x,x + x^x\y,y + y]x\z,z + z\, the corresponding 
set of vertices is denoted by 


N{K) := {x,x + x} X {y,y-Ky} x {z,z-i-2}. 

We refer to the elements of Af := IJ N{K) as nodes. We define fhe active region 

M ;= [rf i.x - r^i] X [[f 1. r - [f i] x [\q-12- [f 1] 

and fhe sef of active nodes Na '■= N 

Definition 4.2 (Skeleton). Given a mesh Q, denote the union of all closed v-orthogonal element faces 
by 5x := {jxeg Sx(^), with 

Hx(^) := {x,x + x} X \y,y + y\ x [z,z + z] 
for any K = [x,x + x]x[y,y + y]x [z, z + z]^§. 

We call 5x the x-orthogonal skeleton. Analogously, we denote the y-orthogonal skeleton by 5y, and 
the z-orthogonal skeleton by E^. 



Definition 4.3 (Global index sets). For any x,y,z € R, we define 

X(y,z) := {x € [0,1] I {x,y,z) e Hx), 

Y{x,z) := |y e [0,T] | {x,y,z) e Hy), 

Z{x,y) := |2 e [0,Z] | {x,y,z) e E^}. 

Note that in an admissible mesh, the entries |0,..., [- k 1 - [+ k ■ ■ ■, 1} are always included 
in X(y, z) (and analogously for Y(x, z) and Z{x, y)). 

Definition 4.4 (Local index vectors). To each active node v = (vi, V2, V3) € Na, we associate a local 
index vector x(v) € which is obtained by taking the unique pi - 1 - 2 consecutive elements in 

X(v2, V 3 ) having vi as their ^^^-th (this is, the middle) entry. We analogously define y(v) e R^2+2 
z(v) e RP3+2_ 




Definition 4.5 (T-spline blending function). We associate to each active node v e Na a trivariate B- 
spline function, referred as T-spline blending function, defined as the product of the B-spline functions 
on the corresponding local index vectors, 

By(x,y,z) := A'xmC-^) • A^y(v)(y) • Nz{v)(z). 


5. Analysis-Suitability 

In this section, we give an abstract definition of Analysis-Suitability. Instead of using T-junction ex¬ 
tensions as in the 2D case, we define perturbed regions fhrough fhe infersecfion of particular T-spline 
supporfs. Analysis-Suifabilify is fhen defined as the absence of intersections between these perturbed 
regions. This idea is comparable to the 2D case, where Analysis-Suitability is defined as fhe absence 
of infersecfions befween T-juncfion exfensions. Subsequenf fo fhese definifions, we prove fhaf all pre¬ 
viously defined admissible meshes are analysis-suifable. 

Definition 5.1 (Perturbed regions). For q,r,s€'K define fhe slices 

SAq) '■= {{x,y,z) e \x = q}, 

Sy{r) := {{x,y,z) e | y ^ r), 

Szis) := {{x,y,z) e 12 = .?). 

Moreover, we denofe by 

A4(<?) := {(vi,V2,V3) € Na I (^Ai.vs) e Hx) 

fhe sef of all nodes of which fhe projecfion on fhe slice S^^iq) lies in some elemenf’s face. Define 
analogously 

Ny{r) := {(vi,V2,V3) e Na I (vi,r,V3) e Hy}, 

A4(5) := {(vi,V2,V3) € Na I {vi,V2,s) € Ez]. 

For any q,r, s eRwe define slice perturbations 

n^{q) := S^{q) n SUppBy n SUppBy, 

veNAq) veNA\NAq) 

'Ryir) := Sy{r) n SUppBy n SUppBy, 

veNy(r) veNA\Ny(r) 

Kzis) := A^zC^) n [J SUppSy n SUppBy. 

veN^is) vsNaWAs) 

The perturbed regions Ry, Rz are defined by 

R. := U := U •= U 

qeR reR isR 

In a uniform mesh, fhe perfurbed regions are empfy. In a non-uniform mesh, fhe perturbed regions are 
a supersef of all hanging nodes and edges (fhis is, all kinds of 3D T-juncfions). See Figure]^ for a 2D 
visualizafion of these definitions. 


9 



Definition 5.2 (Analysis-suitability). A given mesh Q is analysis-suitable if the above-defined per¬ 
turbed regions do not intersect, i.e. if 

= Ky = 0 . 

The set of analysis-suitable meshes is denoted by PEP. 

Remark. When applied in the two-dimensional case, the above definitions may yield perturbed regions 
that are larger than the T-junction extensions from (H [9l (see Fig. [^. However, this occurs only in 
meshes that are not analysis-suitable, and the 2D version of Definition [5]^ is, regarding refinements of 
tensor-product meshes, equivalent to the classical definition of analysis-suitability. 


3^ 



Figure 7: 2D example for the construction of the slice perturbation 7?x(^) in an analysis-suitable mesh. 

The left figure illustrafes fhe consfrucfion of Nxiq) and ifs complemenf Na \ A/x(<?), and fhe 
righf figure shows fhe resulfing slice perfurbafion, which coincides wifh fhe corresponding 
classical T-juncfion exfension. 


Theorem 5.3. AP’'” c for all m> 2. 

Proof. We prove fhe claim by induction over admissible subdivisions. Assume K e Q € A^’'"r\PE'^ and 
lef ^ := subdivf^, A) e AP’™ be an admissible subdivision of We have fo show fhaf ^ e jASP. We 
assume wifhouf loss of generalify fhaf £{K) = 0 mod 3. Hence subdividing K adds m - I faces fo fhe 
mesh, which are x-orfhogonal. SefA [x, x-i-.?]x[y,y-i-y]x[z, z-i-2] andH := {x-\-^x \ j e {1,..., m-1)}, 
then the skeletons of ^ are given by 


Hx = U H X [y,y-i-y] X [z,z-i-z], 5y = 5y, 5^ ^ 5^. 


Let V e Na \ Na be a new active node. Using the local quasi-uniformity from Lemma 


3.5 


it can be 


verified fhaf for all r € 1. such fhaf v € Ny{r) follows Ryir) n supp Bp = 0. Consequenfly, Ky = Ky and 
analogously Rx = Rz- Moreover, R^iq) - Rxiq) for all <7 ^ S. If remains fo characferize 

supp Bv n supp By 

vsAxCf) veNA\NA^) 


Rxi^) - SxiO n 
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Figure 8: 2D example for the construction of the slice perturbation “RxC*?) in a mesh that is not analysis- 
suitable. The left figure illustrates the construction of NJ^q) and its complement Na \ 
and the right figure shows fhe resulfing slice perfurbafion, which is sfricfly larger fhan fhe 
corresponding classical T-juncfion exfension. 


for any ^ € 5. Wifh 


if follows 


NA^) = NA^) ^Na\Na 

and Na \ NA^) = Na \ NAO = Na \ NAO, 


( 2 ) 


- .Sx(^) n \J SUppBy n SUppBy, 

veyVx© veNA\A{^) 

@S.(f)n(lJ SUpp By U SUpp Byj n SUpp By 

veAtxCf) v£A6i\A6i veNaWAO 

= KAO u (SAO n \J suppB{i n suppBy^. 

veKaM^a veNa\NA^) 

'' --’ 

z 


We will prove below fhaf I.r]K.^ - 'LnKy - d). See Figures|^and 10 for an example wifh £(K) = 3 and 
m - 2. Assume for confradicfion fhaf fhere is 5 € K. wifh Kz(s) OE ^ 0. Then fhere exisf v € Nz(s) and 
w € Na \ Nz(s) such fhaf 

Sz(s) nsupp By n supp Bw nE ^ 0. (3) 

Since fhe subdivision of K is admissible, we know fhaf all elemenfs in are af leasf of level {{K)- 

This implies fhaf all fhose elemenfs are of equal or smaller size than K. Denote mid(/f) (cr, v, r) and 

ScU^P’^CA), (4) 

Na \ Na c [cr - cr -i- e] x [v - e, v -i- e] x [t - t -i- e], 


e := It follows 


and with 
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we get more precisely 


S c {^} X [y - £{p 2 + 2), V + s{p 2 + 2)] x [t - £(773 + 2), r + £(p 3 + 2)]. (5) 

The second-order patch := [jK'egp '”iK) consists of elements that may be larger 

in z-direction, but are of same or smaller size than K in ;c- and y-direction. For >v = (wi,W 2 ,>V 3 ), 
Equation Q implies supp and we conclude from Q that 

(>V 1 ,W 2 ) e [^-£(/ 7 i -I- l),^-i-£(/ 7 i -I- 1)] X [y-£( 2 p 2 + 3),v-r £( 2 p 2 + 3)] ( 6 ) 

We assume that there is no element in Q with level higher than ({K) + 1. This is an eligible assump¬ 
tion, since every admissible mesh can be reproduced by a sequence of level-increasing admissible 
subdivisions; see ifTTl Proposition 4.3] for a detailed construction. This assumption implies that the 
z-orthogonal skeleton 5z is a subset of the z-orthogonal skeleton of a uniform {i{K) + l)-leveled mesh, 

5z(@) c Sz(^„|f(^)+i), (7) 

and with = i{K), we have even equality on the patch 

(S) 

using the notation Hz(@P’'”(.K')) := Sz(@) n IJ Since v € we know that Nzi^) ^ 

which means that there are elements in 0 that have z-orthogonal faces at the z-coordinate s, i.e., 
Sz{s) n 5z(^) * 0- With 0 we get Sz{s) n SzC^H^m+i) 0- Since 0u\{(K)+i is a tensor-product 
mesh, its z-orthogonal skeleton consists of global domain slices, which yields Szis) c Ezi0u\e{K)+i)- 
The restriction to the patch 0^’"\K) yields 

Szis) n U 0^'’"iK) C ^ c Ez(0). (9) 

Equation Q implies that Szis) n E 0, and with Q we get that Szis) n IJ 0^’"\K) 0. Hence 

.SzWnU@P’'"(^) 3 Szis)nU^’'”iK) 

= eilpi -I- 3),^ -I- silpi + 3)] X [y - si2p2 - 1 - 3),y - 1 - £(2p2 + 3)] x { 5 }. (10) 

Since w ^ Nzis), we know by definition that (wi, W 2 , 5 ) ^ 5z. Then it follows from 0 that (wi, W 2 , s) ^ 
Szis) n IJ 0^’'^iK), and hence 

(wi,W 2 ) ^ [4^-£(2pi -I- 3),^-i- £(2/7i -r 3)] X [y-£(2p2 + 3),y -I- £(2p2 + 3)] (11) 

in contradiction to @. This proves that *Rz 01, = Q. Similar arguments prove that X D 7?y = 0, which 
concludes the proof. □ 


6. Dual-Compatibility 

This section recalls the concept of Dual-Compatibility, which is a sufficient criterion for linear inde¬ 
pendence of the T-spline functions, based on dual functionals. We follow the ideas of ifTOl for the 
definitions and for the proof of linear independence. In addition, we prove that all analysis-suitable 
(and hence all admissible) meshes are dual-compatible and thereby generalize a 2D result from |5i 
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Figure 9: jz-view on the slice The numbers denote element levels, and the element in the center 

with level 4 is a child of K. The patch is highlighted in blue, and the second-order 

patch is indicated by a thick blue line. 


Proposition 6.1 (Dual functional, lfT3l Theorem 4.4111. Given the local index vector X = (vi,..., Vp+ 2 ), 
there exists an I?'-functional Ax with supp Ax = supp Nx such that for any X - {x\,..., Xp+ 2 ) satisfying 

W xe{xu...,Xp+ 2 } : xi < x < Xp +2 ^ x e {xu... ,Xp+ 2 } 
and d X ^ [x\,... ,Xp+ 2 }'■ x\ < x < Xp +2 => x ^ [x\,... ,Xp+ 2 ], 

follows AxiNx) = Sxx- 

Proof Following |[T3ll . we construct a dual functional on the same local knot vector X which we denote 
by Ax : T^([0,1]) —> R. For details, see llT3l Theorem 4.34, 4.37, and 4.41]. Let yy = cos(^^^7r) for 
j - 0,... ,p + \. Using divided differences, the perfect B-spline of order p - 1 - 1 is defined by 

:= (p -H 1) (-l)P'^n3'o, ■ • • ((Jc - •)+y 

and satisfies (amongst other things) B*^^fx) dv = 1 as depicted in Figure 11 Set 
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Figure 10: 3 ^z-view on the slice Sxi^)- 'Rx is indicated by red areas. Ry is depicted by horizontal red 
lines, R 2 are vertical red lines. At the same time, the squared red area in the center coincides 
with Z. 


Gxix) := 


-L 


2x-Xi-Xp+2 

Xp+2-Xl 


Rp+i(0 dt for xi < X < Xp +2 


and 


<Px(x) = ■1(X-X2)---(X- Xp+l). 


We define the dual functional by 

rxp+2 

Mf)= fDP^\Gx(px)dx forall/eL2([0,l]). 

Jx\ 


(13) 


Note in particular that for all / € L^(R) with f\[xuXp+ 2 ] - 0 follows Axif) - 0. If ( [T^ holds then the 
claim follows by construction, see ifT^ Theorem 4.41]. □ 


We say that two index vectors verifying ( fT^ overlap. In order to define the set of T-spline blending 
functions of which we desire linear independence, we construct local index vectors for each active 
node. 


Definition 6.2. We define the functional Ty by 


Av{Bw) •“ dx(v)(Nx(>v)) ■ dy(v)(Ay(vy)) ' Ax(v){Nz{w)} 

using the one-dimensional functional Ax defined in ([T^. 
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Figure 11: Plot of the perfect B-splines (solid), (dotted), 5*^ (dashed) and the corresponding 
antiderivatives. 


Definition 6.3. We say that a couple of nodes v, w € A/' partially overlap if their index vectors overlap 
in at least two out of three dimensions; this is, if (at least) two of the pairs 

(x(v),x(w)), (y(v),y(w)), (z(v),z(>v)) 

overlap in the sense of Proposition |6.1| 


Definition 6.4. A mesh Q is dual-compatible (DC) if any two active nodes v,w e Na with |supp By n 
supp Bw\ > 0 partially overlap. The set of dual-compatible meshes is denoted by DC^. 


Remark. The above Definition 6.3 fulfills the definition of partial overlap given in ifTOl Def. 7.1], 
which is not equivalent. The definition given in ifTOl is more general, and the corresponding mesh 


classes are nested in the sense DC c ITug). However, we do have equivalence of these definitions in 
the two-dimensional setting. 

The following lemma states that the perturbed regions from Definition |5 .1 1 indicate non-overlapping 
knot vectors, and it is applied in the proof of Theroem| 6 . 6 |below. 


Lemma 6.5. Let q € [0,X] and vi, V 2 € Na- Ifv\ € Nxiq) ^ V 2 and S^iq) n suppByi D suppBv 2 ^ 0. 
then x(vi) and x(v 2 ) do not overlap in the sense of 


This holds analogously for Ny(r), r € [0, F] and Nzis), s € [0, Z]. 


Proof. Let vi = {xi,y\,z\). From vi € Afx(^) and Definition 113 we conclude that {q,y\,z{) e 5x, 
and hence q e X(yi,zi). Let x(vi) = (.vj,..., be the local x-direction knot vector associ¬ 
ated to vi, then supp By, cS^iq) + 0 implies that x\ < q < This and q € X(yi,zi) yield 

q e x(vi). Let V 2 = {x 2 ,y 2 ,Z 2 ). From V 2 ^ Nxiq), we get (q,y 2 ,Z 2 ) ^ ^x, hence q ^ X(y 2 ,Z 2 ), 
and in particular q ^ x(v 2 ). Let x(v 2 ) = (x!^, ..., ) be the local knot vector associated to V 2 , then 

supp Byj n Sxiq) 7 ^ 0 implies that < <7 < Together with x(vi) 3 q ^ x(v 2 ), we see that vi and 

V 2 do not overlap. □ 
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Theorem 6.6. - ITP. 


Proof, “c”. Assume for contradiction a mesh Q which is not DC, hence there exist active nodes 
V, w € Na with |suppBv n suppS^I > 0 that do not overlap in two dimensions, without loss of gener¬ 
ality X and y. We will show that there exist two slice perturbations and 'Ryir) with nonempty 

intersection. We denote v = (vi, V2,V3), w = (wi,W2,W3) and x(v) = {x\,... The elements of 

y(v),x(w),y(w) are denoted analogously. Moreover we define 

Xm := max(xj,x7), xm '■= inin(Xp^^2’-^p,+2) 

:= max(y';,yp, yw := min(y);^^2’3'p2+2) 

Zm := max(zfzi), Zm '■= inin(zJ,^_^2’^p3+2) 

and note that 


SUppB,, n SUpp^H, ^ [Vm,XM] X [yniDM] X [Zm,ZM]- 

Since x(v) andx(>v) do not overlap, there exists q € [xm, xm] with either x(v) 3 q ^ x(w) orx(v) ^ q e x(>v). 
Without loss of generality we assume x(v) 3 q ^ x(w). Since x(v) c X(v 2 , V 3 ), it follows by definition 
that iq,V 2 ,vj) e 5x and hence v e Nyfq). Since q ^ x(w) and hence (< 7 , W 2 ,W 3 ) ^ Hx, it follows that 
w ^ Nxiq)- Then 

n^{q) = S„iq) n suppB,,' n suppB,,' 

v'eN^ig) v'eNA'.N.iq) 

2 Sjiiq) n supp By n supp Bw 
= [q] X [ym,yM] X [Zm,ZM]- 


Analogously, we have 


Py{r) 2 [Xm,VM] X {r} X [Zm,ZM] 


and hence 


P^iq) n 'Ryir) 2 {q} X {r) X [Zm,ZM] * 0, 

which means that the mesh ^ is not analysis-suitable. 

“2”. Assume for contradiction that the mesh is not analysis-suitable, and w.l.o.g. that there is 
V = iq, r, s) € such that R^f^Ry ^ {v) 0. Definition |5 .1 1 implies that there exist vi, V 2 , V 3 , V 4 € Na 

with vi e Nxiq) i V 2 and V 3 € Afy(r) ^ V 4 such that 

V e Sxiq) n Syir) n supp By^ n supp By.^ n supp By^ n supp By^. 

Lemma [O] yields that x(vi) and x(v 2 ) do not overlap, and that y(v 3 ) and y(v 4 ) do not overlap. 

Case 1. If vi e Ny{r) ^ V 2 , or vi ^ Afy(r) 9 V 2 , then vi and V 2 do not partially overlap. 

Case 2. If vi e Afy(r) and V 4 ^ Nxiq), then vi and V 4 do not partially overlap. 

Case 3. If v\ ^ Nyir) and V 3 ^ Nxiq), then vi and V 3 do not partially overlap. 

Case 4. If V 2 e Nyir) and V 4 € Nxiq), then V 2 and V 4 do not partially overlap. 

Case 5. If V 2 4- Nyir) and V 3 € Nxiq), then V 2 and V 3 do not partially overlap. 

In all cases (see the table in Fig. [T^, the mesh is not dual-compatible. This concludes the proof. □ 
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Vi e Afy(r) 

V2 € Nyir) 

V3 € Nyir) 

V 4 e Nyir) 

case(s) 

true 

true 

true 

true 

4 

true 

true 

true 

false 

2 

true 

true 

false 

true 

4 

true 

true 

false 

false 

2 

true 

false 

true 

true 

1,5 

true 

false 

true 

false 

1,2,5 

true 

false 

false 

true 

1 

true 

false 

false 

false 

1,2 

false 

true 

true 

true 

1,4 

false 

true 

true 

false 

1 

false 

true 

false 

true 

1,3,4 

false 

true 

false 

false 

1,3 

false 

false 

true 

true 

5 

false 

false 

true 

false 

5 

false 

false 

false 

true 

3 

false 

false 

false 

false 

3 


Figure 12: The five cases considered in the proof of Theorem 6.6 cover all possible configurations. 


Theorem 6.7. Let Q be a DC T-mesh. Then the set of functionals {Ay \ v e Na) D a set of dual 
functionals for the set [By \ v € Na)- 

The proof below follows the ideas of 191 Proposition 5.1] and ifTOl Proposition 7.3]. 

Proof Let v,w € Na- We need to show that 

AyiByy) = 6yyy, (14) 


with 5 representing the Kronecker symbol. 

If supp By and supp By, are disjoint (or have an intersection of empty interior), then at least one of 
the pairs 


(supp(A/x(v)), supp(A/x(w))), (supp(A/y(v)), supp(77y(^v))), (supp(77z(v)), supp(A^z(n,))) 

has an intersection with empty interior. Assume w.l.o.g. that |supp(Ax(v)) supp(Ax(w))| - 0, then 

Av{By,} = Ax(y)(-lVx(w)) ■dy(v)(A^y(w)) ■ '^z(v)(l^z(w)) ~ 0. 

0 

Assume that supp By and supp By, have an intersection with nonempty interior. Since the mesh ^ is 
DC, the two nodes overlap in at least two dimensions. Without loss of generality we may assume the 
index vectors (x(v),x(w)) and (y(v),y(w)) overlap. Proposition |6.1| yields 

Ax(v)(Nx(n,)) — and Ay(y)(Ny(y,)) — 

^V2W2' 

The above identities immediately prove ( fTT] ) if vi wi or V 2 t W 2 - If on the contrary, vi - wi 
and V 2 t W 2 , then v and w are aligned in z-direction, this is, z(v) and z(>v) are both vectors of p + 2 
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consecutive indices from the same index set Z(vi, V 2 ) = Z(wi, W 2 ). Hence v and w must overlap also in 
z-direction. Again, Proposition|6.1|yields 


'^Z(v)(Az(h>)) — (5^3 VP 3 ? 


which concludes the proof. □ 

Corollary 6.8 t llTOl Proposition 7.4]). Let Q be a DC T-mesh. Then the set {B^ \ v € Na) is linear 
independent. 


Proof. Assume 


^ CyB,, = 0 
veNa 


for some coefficients {cvlveAOi £ Then, for any w e Na, applying A„ to the sum, using linearity and 
Theorem |6.7[ we get 

Cyy — ^ ^ CyBy^ — 0 . □ 

veNa 


7. Linear Complexity 


This section is devoted to a complexity estimate in the style of a famous estimate for the Newest Vertex 
Bisection on triangular meshes given by Binev, Dahmen and DeVore lfT4l and, in an alternative version, 
by Stevenson |[T5]| . Linear Complexity of the refinement procedure is an inevitable criterion for optimal 
convergence rates in the Adaptive Finite Element Method (see e.g. ifTdlfTSlfT^ and ifTTl Conclusions]). 
The estimate and its proof follow our own work IfTTl[TSll . which we generalize now to three dimensions 
and m-graded refinement. The estimate reads as follows. 

Theorem 7.1. Any sequence of admissible meshes Qq, Qi,...,Qj with 

Qj = refP’™!^,-!, 7Vf,_i), 7Vf,_i c for j € [\,... ,J} 


satisfies 


with Co m - 


7-1 


„l/3 


P-m - 1 - 


-1/3 


\Qj\Qo\ ^ Cp^m 'y, \Mj\ , 

7=0 

{Ad\ + 1) ( 4 r /2 +m^i^) {Adj, + ni^i^) and di,d 2 ,dj, from Lemma 


7.2 


below. 


Lemma 7.2. Given M. Q Q ^ 
i{K) < {{K') + 1 and 


and K e refP’'”(^, vVl) \ Q, there exists K' ^ M. such that 


with “< ” understood componentwise and constants 


d\ 

d2 - ^ 


d^ 


\-m 

1/3 


{pi + 
{P 2 + 
(P3 + 


2 rrp- ' 

2 rrP- ' 

o “I" ..D 


The proof is given in Appendix |A.3| 
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Proof of Theorem \7.1\ 

(1) For A' € IJ AP’™ and A € At := Mq U • • • U Aty-i, define A{K, K) by 

i if + f Dist(A, k) < , d 2 , df), 

A{K,k) := 

10 otherwise. 

(2) Main idea of the proof 

\Gj\Qo\= 2] 1 f 2] 2]d(A,A) 

Kegj\go KeQpQo KeM 

|4l ^ 

s Yj Cp,m = Cp,m Y 
KeM ;=0 


(3) Each K e Qj\Qq satisfies 


Y 


> 1 . 


KeM 


Consider K € Qj\Gq. Set j\ < J such that K € Gj^+\\Gji ■ Lemma 7.2 states the existence of K\ € Mji 
with Dist(A, Ai) < m--^^^'>>^{di,d 2 ,d 3 ) and {(K) < e{K\) + 1. Hence A{K,Ki) = > 0. The 

> 72 > 73 > ■ ■ ■ and A 2 , A 3 ,... with A-i e Qj.+i \ §j^ and 


repeated use of Lemma |7i2] yields ji 
Ki e Aly, such that 


Dist(A,_i, A;) < m-^^^‘-^'^l^{dud 2 ,d 3 ) and AA;-i) < ^(A,) + 1. 


(15) 


We repeat applying Lemma 7.2 as A{K,Ki) > 0 and 7’(A,) > 0, and we stop at the first index L with 
A{K, Kl) - 0 or 7 ’(Al) = 0. If €{Kl) - 0 and A{K, Kl) > 0, then 


Y KK,k) > A{K,Kl) 


= m 


it{K)-e(K[,))i^ > ^ 1/3 


KeM 


If A{K, Kl) = 0 because £{K) > {{Kf) + 1, then ( [T5] ) yields {{Ki-i) < i{Ki) + 1 < i{K) and hence 

Y ^iK,k) > d(A,AL_i) - ^(W-d/fL-i))/3 > ^1/3 


KeM 


If A{K, Kl) = 0 because Dist(A, Kf) > Im ^^^^^^{di,d 2 , d^), then a triangle inequality shows 

L-l 

2;„-W/3(Ji J2, J3) < Dist(A, Ai) + 2 ] Dist(A, A+i) 

1=1 

< d2, d^) + ^ d2, d^). 


i=l 


and hence m' 


-£(K)I3 


L-l 

<Ym 

i=\ 


-q/r,)/3 


L-l 


. The proof is concluded with 


L-l 


1 < ^ < Y A{Kk). 


1=1 


i=i 


KeM 
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(4) For all j e {0,. .. ,J — 1) and K € Mj holds 

Yj ^ + 1) (4^2 + {4d3 + m2/3) ^ Cp,^ . 

Kegj\go 

This is shown as follows. By definition of A, we have 
2] AiK,K)< Y 

Kegj\go Ke\JAf’"'\go 

e(K)+l 

= Y # {^ e U AP’™ I e{K) = j and Dist{K, K) < 2m-j^\di ,d 2 ,d 3 )}. 

A* " B ' 

Since we know by definition of the level that £{K) = j implies |.K'| = m~j, we know that | 1 J B| is an 
upper bound of #B. The cuboidal set IJ B is the union of all admissible elements of level j having their 
midpoints inside a cuboid of size 

4m~^^^d\ X 4m~^^^d2 x 4m~^^^d3. 


An admissible element of level 7 is not bigger than x x Together, we have 

|lj B| < m~-^ {4d\ + 1) (4^2 + {4d3 + 

and hence #B < (4di + 1) {4d2 + (4^3 + An index substitution k \= \ - j + £{K) proves the 

claim with 


i(K)+\ e(K) 

7=1 t :=0 




< m 


1/3 


Z' 

t :=0 


-k/3 


m‘/3 


□ 


An experiment on C 


p,m 


1 smaller ratios of 


131. Starting from 


The constant Cp m arising from this theory is very large, however we observed muc 
refined and marked elements in the experiment (in all cases less than , see Figure 
a 5 X 5 X 5 mesh, we applied the refinement algorithm with only one corner element marked, always 
sticking to the same corner. This is realistic when resolving a singularity of the solution of a discretized 
PDF. The advantage of greater grading parameters could not be seen in random refinement all over the 
domain. 


8. Conclusions & Outlook 

We have generalized the concept of Analysis-Suitability to three-dimensional meshes that originate 
from tensor-product initial meshes, and proved that it guarantees linear independence of the T-spline 
blending functions. We introduced a local refinement algorithm with adjustable mesh grading, and 
proved that it has linear complexity in the sense that the overhead for preserving Analysis-Suitability 
is essentially bounded by the number of marked elements. We expect that these results also generalize 
to even-degree and mixed-degree T-splines. In order to achieve this, a universal definition of anchor 
elements is needed, based on the techniques from ifT^ . 

Open questions that have not been investigated in this paper address the overlay (this is, the coarsest 
common refinement of two meshes), the nesting behavior of the T-spline spaces, and more general 
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grading parameter m 



grading parameter m 


Figure 13: The complexity constant Cp m in theory (left) and experiment (right). The values of C'„ were 
taken from an experiment illustrated in Figure [T4| 


meshes. As in our preliminary work ifTTI . we expect that the overlay has a bounded cardinality in terms 
of the two overlaid meshes, and that it is also an admissible mesh. Nestedness of T-spline spaces is not 
evident in general EOl . but we expect nestedness for the meshes generated by the proposed refinement 
algorithm. A first step in this issue will be a characterization of three-dimensional meshes that induce 
nested T-spline spaces. A generalization of this paper to a more general class of meshes will most likely 
require a manifold representation of the mesh, and use recent results on Dual-Compatibility in spline 
manifolds ll^ . 


References 

[1] T. Sederberg, J. Zheng, A. Bakenov, and A. Nasri, T-Splines and T-NURCCs, ACM Trans. Graph. 
22 (2003), no. 3, 

[2] D. R. Forsey and R. H. Bartels, Hierarchical B-spline refinement, Comput. Graphics 22 (1988), 
205-212. 

[3] G. Kuru, C. Verhoosel, K. van der Zee, and E. van Brummelen, Goal-adaptive isogeometric 
analysis with hierarchical splines, Comput. Methods Appl. Mech. Engrg. 270 (2014), 270 - 292. 

[4] C. Giannelli, B. Juttler, and H. Speleers, THB-splines: the truncated basis for hierarchical 
splines, Comput. Aided Geom. Design 29 (2012), 485^98. 

[5] T. Dokken, T. Eyche, and K. Pettersen, Polynomial splines over locally refined box-partitions, 
Comput. Aided Geom. Design 30 (2013), no. 3, 331 - 356. 

[6] E. J. Evans, M. A. Scott, X. Ei, and D. C. Thomas, Hierarchical T-splines: Analysis-suitability, 
Bezier extraction, and application as an adaptive basis for isogeometric analysis, Comput. Meth¬ 
ods Appl. Mech. Engrg. 284 (2015), 1-20. 


21 





Figure 14: Estimation of the experimental constants C', for m = 2,..., 5. 


[7] A. Buffa, D. Cho, and G. Sangalli, Linear independence of the T-spline blending functions associ¬ 
ated with some particular T-meshes, Comput. Methods Appl. Mech. Engrg. 199 (2010), no. 2324, 
1437 _1445. 

[8] X. Ei, J. Zheng, T. Sederberg, T. Hughes, and M. Scott, On Linear Independence of T-spline 
Blending Functions, Comput. Aided Geom. Des. 29 (2012), no. 1, 63-76. 

[9] E. B. da Veiga, A. Buffa, D. Cho, and G. Sangalli, Analysis-Suitable T-splines are Dual- 
Compatible, Comput. Methods Appl. Mech. Engrg. 249-252 (2012), 42-51, Higher Order Einite 
Element and Isogeometric Methods. 

[10] E. B. da Veiga, A. Buffa, G. Sangalli, and R. Vazquez, Mathematical analysis of variational 
isogeometric methods, Acta Numerica 23 (2014), 157-287. 

[11] R Morgenstem and D. Peterseim, Analysis-suitable adaptive T-niesh refinement with linear com¬ 
plexity, Comput. Aided Geom. Design 34 (2015), 50-66. 

[12] M. Scott, X. Ei, T. Sederberg, and T. Hughes, Local refinement of analysis-suitable T-splines, 
Comput. Methods Appl. Mech. Engrg. 213216 (2012), 206-222. 

[13] E. Schumaker, Spline Functions: Basic Theory, 3 ed., Cambridge Mathematical Eibrary, Cam¬ 
bridge Univ. Press, Cambridge, 2007. 

[14] P. Binev, W. Dahmen, and R. DeVore, Adaptive Finite Element Methods with convergence rates, 
Numer. Math. 97 (2004), no. 2, 219-268. 

[15] R. Stevenson, Optimality of a standard adaptive finite element method. Pound. Comput. Math. 7 
(2007), no. 2, 245-269. 

[16] C. Carstensen, M. Peischl, M. Page, and D. Praetorius, Axioms of adaptivity, Comput. Math. Appl. 
67 (2014), no. 6, 1195-1253. 


22 















[17] A. BufFa and C. Giannelli, Adaptive isogeometric methods with hierarchical splines: error esti¬ 
mator and convergence, ArXiv e-prints (2015). 

[18] A. BufFa, C. Giannelli, P. Morgenstem, and D. Peterseim, Complexity of hierarchical refinement 
for a class of admissible mesh configurations, Computer Aided Geometric Design (2016), -, In 
Press. 

[19] L. B. da Veiga, A. BufFa, G. Sangalli, and R. Vazquez, Analysis-suitable T-splines of arbitrary 
degree: definition, linear independence and approximation properties. Math. Models Methods 
Appl. Sci. 23 (2013), no. 11, 1979-2003. 

[20] X. Li and M. A. Scott, Analysis-suitable T-splines: Characterization, refineability, and approxi¬ 
mation, Math. Models Methods Appl. Sci. 24 (2014), no. 06, 1141-1164. 

[21] G. Sangalli, T. Takacs, and R. Vazquez, Unstructured spline spaces for isogeometric analysis 
based on spline manifolds, ArXiv e-prints (2015). 


A. Minor proofs 

A.1. Proof of Lemma |33l 

If K - K, the claim is trivially Fulfilled. IF otherwise K K,\\/e consider the Following two cases. 

Case 1. Assume that £{K) - £{k) 1. Since K = [x,x x\y:\y,y + y^x\z,z + z\ is the result oF 

successive subdivisions oF a unit cube, it holds that 




size(^(X)) := {x,y,z) - ■ 


,„-(r(^)-2)/3(l 1 1) 


iF €{K) - 0 mod 3, 
iF €{K) - 1 mod 3, 
iF i{K) - 2 mod 3. 


Since K results From the subdivision oF K, we also have that 


'(m-(^m+«/3,o,0) 


Dist(X,X) - - 


(0,m-W«+5)/3^0) 


(0,0,m-(W+4)/3) 


iF i{k) = 0 mod 3, 
iF i{k) - 1 mod 3, 
iF i{k) - 2 mod 3. 


Recall that 




DP’'"(k) := 


We rewrite ([TT]) in the Form 


(Pl + l,P2 + |,P3 + i) 


m 


Pi + i,P 3 + i) 


ft-l)/3 / Pi+3/2 , 3 

m 

(A:-2)/3(- Pi+3/2 P2+3/2 , 3^, 

V. m ’ m ^ 2 ' 


ifk = 0 mod 3, 
ifk = \ mod 3, 
ifk = 2 mod 3. 


Dist(X,X) 


'(0,0,m-W^)+3)/3) 

(m-(^m+ 5 )/ 3 ,o, 0 ) 

(0,m-W^)+4)/3^0) 


iF i{K) - 0 mod 3, 
iF i{K) - 1 mod 3, 
iF i{K) - 2 mod 3 


(16) 


(17) 


(18) 
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and observe that DP’'”(^(^)) + Dist(^, K) < - 1) = DP’'”(^(^)). The case 1 is concluded with 


= {v e n I Dist(^,x) < DP’^C^C^))) 
c {;c e a I Dist(^, x) < DP’'”(^(^)) + Dist{K, K)} 
c U^’”XK), 


and consequently gP’'”{K) c gv^^{k). 

Case 2. Consider K <z K with {{K) > £{K) + 1, then there is a sequence 

K = KocKiC---cKj = k 

such that Kj-i € child(^y) for j - 1,..., L. Case 1 yields 

@P’'”(/:) c c... c @P’'”(^). □ 

A.2. Proof of Lemma l33] 

For £{K) = 0, the assertion is always true. For €{K) > 0, consider the parent K of K (i.e., the 
unique element e IJ AP’'” with K e child(A)). Since g is admissible, there are admissible meshes 
g^,... ,gj = g and some y € {0,...,7 - 1) such that K € gj+x - subdiv(^y, {A}). The admissibility 
gj+x e AP’™ implies that any K' e g^j'”{k) satisfies €{K') > €{k) = £{K) - 1. Since levels do not 
decrease during refinement, we get 


i{K) - 1 < min7(@P’'”(A)) < mm£{g^'^{k)) 


(19) 



□ 


A.3. Proof of Lemma |7^ 


The coefficient DP’'”(A:) from Definition |2.4| is bounded by 

DP’'”(/:) < {jjx + \, + |), for all A: € N. 


( 20 ) 


p 

Recall sizefAc) from (fT^ and note that it is decreasing and bounded by 


size(A:) < m (l,m'^^,m^^^). 


( 21 ) 


Hence for A € ^ € AP’™ and k' e g^''’\k), there is v e A' n f7P’™(A) and hence 

Dist(A, k') < Dist(A, x) + Dist(A', x) 



s 


< m (p + s). 


( 22 ) 
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The existence of ^ e refP’'”(^, Al) \ @ means that Algorithm |2.9| subdivides K' - Kj,Kj-\,... ,Kq 
such that Kj-\ e and C{Kj-i) < €{Kj) for 7 = 7,, 1, having K' € M and K € child(Ao), 

with ‘child’ from Definition 2.6 Lemma [33] yields i{Kj-\) - i{Kj) - 1 for 7 = 7,..., 1, which yields 
the estimate 


J 1I22I ^ 

Dist(r,Ao) < 2]Dist(A7,A7_i) T + 


7=1 

j 


,/=i 


= ^ m (p + s) <m (p + s) ^ rn 


■7/3 


7=1 


7=1 


m 


-i/3-mo)P 


l — m 


(P + s) 


m 


<K)/3 


I — m 


(P + s). 


From ([18]) we get 


Dist(Ao,A) < 


-(t(K)+5)l3 ^-(W+4)/3 ^-{e(K)+3)l3 

This and a triangle inequality conclude the proof. 
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